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Abstract. A short review of Algebraic Geometry tools for the decomposition of tensors and polynomials is given 
from the point of view of applications to quantum and atomic physics. Examples of application to assemblies of indis- 
tinguishable two-level bosonic atoms are discussed using modern formulations of the classical Sylvester's algorithm 
for the decomposition of homogeneous polynomials in two variables. In particular, the symmetric rank and symmetric 
border rank of spin squeezed states is calculated as well as their Schrbdinger-cat-like decomposition as the sum of 
macroscopically different coherent spin states; Fock states provide an example of states for which the symmetric rank 
and the symmetric border rank are different. 

PACS. 03.65.Ud Quantum entanglement - 02. lO.-v Algebraic geometry - 42.50. Dv Squeezed states 



1 Introduction 

The mathematical research on tensor decomposition is gener- 
ally considered to have started out in the late XIX century with 
the pioneering work of J. J. Sylvester on the decomposition of a 
homogeneous polynomial in two variables as the sum of pow- 
ers of linear forms. Still now, this kind of questions are attract- 
ing a strong interest from the algebraic geometry community 
and they have been related to the problem of decomposing a 
multi-dimensional, multi-index tensor as the sum of simpler, 
rank-1 objects; the minimum number of elements of such a de- 
composition is called the rank of the tensor. Given the isomor- 
phism between homogeneous polynomials and totally symmet- 
ric tensors, Sylvester's result is a first example in this program. 
In the last years, the algebraic geometry tools for tensor decom- 
position are starting being applied in a number of application 
fields, ranging from Signal Processing CEE], Telecommuni- 
cations J4][5] to Phylogenetics 07 8 1, Chemometrics J9] and 
Algebraic Statistics II10II1 11 . Geometric complexity theory lfT2l 
[T3lfT4l . Arithmetic complexity fll5lll6l[T7l[lin . Data Analysis 
(e.g. Independent Component Analysis [19 20,27] ). 

In physics, the Hilbert space of a composite quantum sys- 
tem is the tensor product of the Hilbert spaces of the constituent 
sub-systems, e.g. the individual particles forming a quantum 
gas or the entangled, yet physically separated sub-systems of 
quantum information processing devices. Also in this context, 
it is often of interest to expand a given quantum state as the 
sum of factorizable states, i.e. tensor products of pure states of 
each sub-system. The minimum number of factorizable states 
that have to be summed to recover a given state is called the 



Schmidt rank of the state and is often mentioned as a possible 
measure of entanglement (22,23|. In many cases, the decom- 
position of the state has itself a physical interest, as it highlights 
an underlying Schrodinger cat structure as a linear superposi- 
tion of macroscopically different states. 

In the present paper, we investigate the application of a al- 
gebraic geometry techniques to the calculation of the rank and 
of the decomposition of quantum states of actual physical inter- 
est of an assembly of N indistinguishable two-level atoms, in 
particular spin squeezed states and Fock states. On one hand, 
two level atoms are amenable to a description in terms of bi- 
variate homogeneous polynomials that can be tackled using 
Sylvester's algorithm and its modern extensions 0241125 ll26l . 
On the other hand, spin squeezed states have been proposed 
as a practical tool to overcome the shot noise limit in preci- 
sion interferometrical measurements [27 ,28,291: experimental 
demonstrations of this idea have been reported using either the 
internal atomic level structure in a cloud of thermal BUI or 
Bose-condensed RD atoms or the atomic center-of-mass mo- 
tion in a double well geometry l32l . 

The article is organized as follows: the relation between 
quantum states, tensors and polynomials is reviewed in Secf2] 
A general introduction to the concept of symmetric rank and 
symmetric border rank is given in SecfJJ The basic definitions 
are given in Sec l3.ll together with a physical interpretation 
in the second quantization formalism: for symmetric (skew- 
symmetric) tensors describing a Bose (Fermi) system, the de- 
composition of the physical state is made over Hartree (Hartree- 
Fock) states. A modern formulation of Sylvester's algorithm of 
ll25l is presented in Sec l3.2l 
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Sec|4] is devoted to the application of these concepts to 
specific physical states, namely spin squeezed states and Fock 
states. The calculation of the symmetric rank and border rank 
of spin squeezed states is presented in Sec l4.lt the symmetric 
rank and symmetric border rank turn out to be always equal 
for this class of states, but interesting features are found when- 
ever the squeezing parameter is a rational multiple of 27r, fi — 
2Trp/q. For these values, the symmetric rank is lower than its 
generic value; as originally predicted in 11331 . the state can be 
decomposed as the sum of q macroscopically different coherent 
spin states pointing on equispaced directions along the horizon- 
tal plane orthogonal to the quantization axis. Sec l4.2l is devoted 
to Fock states with a well defined occupation number of the 
two states, for which the symmetric rank and symmetric bor- 
der rank turn out to be different. Explicit decompositions are 
given as an exact sum or as the limit of sums of coherent spin 
states aligned on the horizontal plane. 

Conclusions and perspectives are summarized in Sec. [5] 



It is immediate to see that this isomorphism can be extended by 
linearity to the whole space of polynomials and to the whole 
many -body Hilbert space H s . With the proper normalization 
factor, the isomorphism between polynomials and physical many- 
particle states induced in second quantization by (0 is fully 
equivalent to the standard one of the mathematical literature 
summarized by (Q3 and (0. 

An analogous isomorphism between physical states and Grass- 
mann-generalized polynomials in anti-commuting variables can 
be defined for systems of fermionic (i.e. with half-integer spin) 
particles for which the many-body wavefunction is totally an- 
tisymmetric under exchange. Such states correspond to skew- 
symmetric tensors and, within second quantization, can be writ- 
ten in the form (0 in terms of creation operators a\ which now 
anti-commute. 



3 Rank and Border rank 



2 Quantum states, tensors and polynomials 

We consider a system of N quantum particles. The Hilbert 
space Hi for each individual atom is assumed to have dimen- 
sion n. The Hilbert space of the A^-body system is obtained 
as the tensor product of N copies of the n-dimensional sin- 
gle particle Hilbert space Hi- In the case of indistinguishable 
bosonic (i.e. with integer spin) particles, the indistinguishabil- 
ity postulate of quantum mechanics imposes that only the to- 
tally symmetric states under particle exchange are physically 
relevant, which amounts to restricting our attention to the sub- 
space H s = S N Hi C Hi of completely symmetric ten- 
sors. 

Using the {ei, . . . , e„} orthonormal basis of Hi, a sym- 
metric tensor t can be written as 

1 < i je{i...N}<™ 

where the coefficients tensor a (i.e. the many-body wavefunc- 
tion in the {ei, . . . , e n } basis) is symmetric under exchange, 

a...^...^,... = a... i . l ...i j ... lor all / ../ , 1 n. 

The standard isomorphism (i.e. biunivocal linear map) be- 
tween symmetric tensors and homogeneous polynomials asso- 
ciates the tensor t E S N H\ defined in (Q} to the following 
homogeneous polynomial p{x\ , . . . , x n ) of degree N, 

p(xi, . . . ,x N ) = 2J ai u ,„ i i N x il ...Xi N . (2) 

1 <Se{i...JV}<n 

A related isomorphism between physical states and polyno- 
mials naturally arises within the second quantization descrip- 
tion of many-body states in terms of the usual creation opera- 
tors a\ for an atom in the i state. Each monomial x™ 1 ■ • • x^f of 
degree A?" is associated to a (unnormalized) physical state with 
a well defined number of particles n, in each single particle 
state i = 1, . . . , n and total number of particles ^2i=i n i = 

x?-x»» (otr-^ot^lvac). (3) 



3.1 Definition and physical interpretation 

For a given symmetric tensor t € S N H\, its symmetric rank 
sr(t) is defined as the minimum integer r for which there ex- 
ists a set of (not necessarily normalized nor orthogonal) vectors 
vi, . . . , v r 6 H\ such that 

*=i>r- (4) 

i=l 

The symmetric border rank sbr(t) of a symmetric tensor t € 
S N Hi is defined as the minimum integer s for which there ex- 
ists a sequence of symmetric tensors ti E S N Hi of symmetric 
rank s tending to t, lirrii_ i . 00 ti = t. Obviously, one has 

sbr(t) < sr(t). (5) 

From a physical standpoint, the symmetric rank sr has a simple 
interpretation in terms of the decomposition of a bosonic many- 
particle state as the sum of Hartree states where all A^ atoms 
share the same single-particle state Vj. The symmetric rank sr 
of a state \ip) is then the minimum number of terms that have to 
appear in the sum in order to exactly reconstruct the state \ip). 

This decomposition can be straightforwardly reformulated 
within second quantization in terms of the associated polyno- 
mial. Given a homogeneous polynomial p(xi, . . . , x n ) of de- 
gree N in n variables, the symmetric rank sr(p) is the mini- 
mum integer r for which there exists a decomposition as the 
sum of r N-th powers of linear forms, 

/ \ N 

r I ri \ 

p{xi,...,x n ) I ^2^i,jXj\ (6) 

»=i \i=i / 

with 4>i,j <G C. 

Identifying again the variables xi,. . . ,x n of the polyno- 
mial with the creation operators oj of the atom in the i = 
1, . . . , n single-particle state, the above definition of symmet- 
ric rank in terms of Hartree states is immediately recovered: it 
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is the minimum integer r for which 




(7) 



The complex coefficients <J>; j of the i-th linear form on the j-th 
variable of the polynomial have the physical interpretation of 
the j-state component of the i-th Hartree wavefunction. 

Analogously the symmetric border rank sbr(p) of a homo- 
geneous polynomial p(xi , . . . , x n ) of degree JV in n variables 
is defined as the minimum integer s for which there exists a se- 
quence of homogeneous polynomials pi (xi , . . . , x n ) of degree 
JV in n variables of symmetric rank s tending to p(xi , . . . , x n ), 
lim i _ ) . 00 p i (a;i, ...,x n )= p(xi, . . .,x n ). 

These definitions can be generalized to not necessarily sym- 
metric tensors, e.g. the ones describing the state of distinguish- 
able quantum particles. For a given tensor t G Hf N , its rank 
r(t) is defined as the minimum integer r for which there ex- 
ists a family of vectors G Hi with j = 1, . . . , JV and 
i = 1, . . . , r, such that 



(i) 



AN) 



(8) 



The border rank br(t) of a tensor t G Hf N is defined as the 
minimum integer s for which there exists a sequence of rank- 
s tensors ti G Hf N such that lim^oo ti = t. Note that in 
the absence of the symmetry requirement, the single-particle 

wavefunctions v^ 1 ' , . . . , of the different N particles are 
no longer necessarily equal. The so-called Comon's conjec- 
ture Oil states that if t e S N Ui, then sr{t) = r(t). Such 
a conjecture is proved to be true only in few cases l34l ?l. An 
analogous conjecture can be formulated in terms of the border 
rank and the symmetric border rank (see e.g. (35)). 

Analogous concepts can be introduced for skew-symmetric 
tensors t G /\ Hi- In this case, the skew-symmetric rank 
ar{t) is defined as the minimum integer r for which there ex- 



ists a family of orbitals v 
i = 1, . . . , r, such that 



G Hi with j = 1,. . ., JV and 



vP A 



A v. 



(JV) 



(9) 



The skew-symmetric border rank abr(t) of a skew-symmetric 
tensor t € /\ Hi is defined as the minimum integer s for 
which there exists a sequence of skew-symmetric rank s skew- 
symmetric tensors ti G f\ N Hi such that lim^^ ti = t. 

In physical terms, the decomposition (O corresponds to ex- 
panding the many-particle state of a Fermi system as the sum 
of Hartree-Fock states, 



E 




k=l 



vac) 



(10) 



1 It is interesting to note that the rank of a general skew-symmetric 
t G /\ 2 Hi is equal to half the rank of the skew-symmetric matrix T 
that is obtained by writing t in an arbitrary basis of Hi 1361 . 



where, <fr^ is the j-state component of the fc-th orbital k = 
1, . . . , N) of the i-th Hartree-Fock state (i = 1, . . . , r) and 
fj are the (fermionic) creation operators for a particle in the 
j single-particle state; as usual for Fermi systems, creation op- 
erators satisfy the anti-commutation rules of Grassmann vari- 
ables, 

(ID 



/7/j-/j// = o 



3.2 Algorithm 



Let now Hi be a 2-dimensional vector space defined over the 
complex numbers. The algorithm to compute the symmetric 
rank of a homogeneous polynomial p in two variables (i.e. the 
symmetric rank of a symmetric tensor t G S N Hi) is classically 
attributed to J. J. Sylvester. A more modern proof of it is pre- 
sented in l24l . We actually refer to a simplified version of the 
same available in l25l . From now on we will interchangeable 
use the polynomial and symmetric tensor language. 
Before going into the details of the algorithm we need to de- 
fine the so-called Catalecticant matrices I371 I38II Cjy-r.rip) 
(Cjv-r,r(t)) associated to a given bivariate polynomial p(x, y) 
(to a given tensor t G S Hi), 



p(x,y) 



N 

1 a lX l y N -\ 



i=0 



(12) 



The ?--th Catalecticant matrix Cat-^Ij?) = (ckj) associated to 
pis the (JV— r+l)x (r + 1) matrix with entries: c^j = ctk+j-2 
with k = 1, . . . , N — r + 1 and j = 1, . . . , r + 1. 

It is classically known l38l that, for a given homogeneous 
bivariate polynomial p(x, y) of degree JV, if r is the minimum 
integer for which all the minors of Cjy-r rip) of order (r + 1) 
vanish then sbr(p) = r. Moreover in fl24l and in l25l it is 
shown that if sbr(p) — r, then sr(p) can only be equal either 
to r or to JV — r + 1. The way to distinguish among these last 
two possibilities is to verify if a certain polynomial of degree r 
constructed via elements in the kernel of CN- r .r{p) (see pre- 
cisely below in the algorithm) has or not distinct roots. So the 
Algorithm that we are going to describe works as follows: for a 
given p(x, y) find firstly its symmetric border rank, then com- 
pute its symmetric rank, and, finally, give the decomposition in 
the form ©. 

It is worth to remark here that, over the complex numbers, 
the symmetric border rank of any homogeneous p(x, y) is al- 
ways bound from above by the symmetric rank of a generic 
bivariate polynomial, that is, 



?br(p) < 



JV + 1 



(13) 



Moreover, since if sbr(p) = r then sr(p) G {r, JV — r + 2}, 
we easily deduce that if p(x, y) is such that its symmetric bor- 
der rank sbr(p) is strictly less than both its symmetric rank 
sr(p) and the generic symmetric rank [^-^1, then the sym- 
metric rank of p(x, y) is strictly bigger then the generic value 

We are now ready to state the Algorithm that for a given 
homogeneous bivariate polynomial of degree JV, calculates its 
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symmetric border rank, its symmetric rank and its decomposi- 
tion as the sum of r TV-th powers of linear forms as in ©. 

Algorithm 

Input: The projective class of a degree TV homogeneous poly- 
nomial p(x, y) in two variables. 

Output 1: sbr(p). 
Output 2: sr(p). 

Output 3: a decomposition of p as in ©. 



1. Initialize r = 0; 

2. Increment r <— r + 1; 

3. Compute CN- r ,r(pYs (r + 1) x (r + l)-minors; 

if they are not all equal to zero then go to step [2] (note that 

this can only happen if r < ); 

else, 

Output 1: sbr(p) = r. Go to step [4] 

4. Choose a solution (uq, . . . ,u r ) of the system 



C N - r ,r(p) ■ (""0, • • • , U r f = 0. 



(14) 



If the polynomial 



g(io,*x)=«o*o + «i*o 1 *x + -"+*tr*i (15) 
has distinct roots, 

r 

q(t ,ti) = /xJJ^fcto - aj.ii), (16) 
fc=i 

then 

Output 2: sr(p) = r and go to step|6l 
otherwise 

Output 2: sr(p) = TV — r + 2 and go to step [5] 
5. Choose a solution (uq, . . . , Tt sr ( p )) of the system 

CiV-sr(t),sr(t)(p) ' ( u 0, • • • > "sr(p))* = 0. (17) 

The polynomial 

/, , \ — ,sr(p) — ,sr(v) — l, . . — j.srlp) 

q(t ,ti) = U t Q + Ul t Q w t\ H h U sr (p)*l 

(18) 



has distinct roots, 



sr(p) 



(19) 



fe=i 



Go to step [6] 

6. Solve the linear system in the variables Ai, . . . , A sr ( p ) given 
by 



sr(p) 



N 



(20) 



k=l 



For such A/s, the decomposition ( f20l > is our Output 3. 



Observe that, for the polynomials p(x, y) for which sbr(p) < 
sr(p), the polynomial q(to, ii) defined in (TT~5T > has not distinct 
roots, i.e. there exist integers si, . . . , s r i with r' < sr(p) and 

Sfe=i s fc = sbr(p) < sr(p) such that the polynomial q(to, t\) 
obtained as in step|4]of our algorithm can be written as follows, 



q{t ,ti) = jj, Jj09fcXo - <5fe2:i) s 



(21) 



fe=i 



We like to stress here that those roots have an algebraic mean- 
ing even if they are not distinct (see l25ll ). In fact, in this case, 
there exist Ai, . . . , \ r > G C and polynomials in two variables 
Fi, . . . , F r i of degrees (si — 1), ... , (s r i — 1) respectively such 
that 



ax 



■Pv) 



N-(s k -l) 



F k . 



(22) 



k=l 



Before moving to the physical applications, we wish to stress 
that in the present section we have chosen to restrict our atten- 
tion to the algorithm that compute the symmetric rank for a 
polynomial in two variables because this is the relevant case 
for the applications that we are going to discuss in the next sec- 
tion. A more general algorithm that computes the symmetric 
rank of a homogeneous polynomials in n > 2 variables can be 
found in ll26ll . 



4 Examples of application 

We now proceed to illustrate the application of the general con- 
cepts discussed in the previous sections to specific examples of 
current interest in the field of atomic physics. In particular, we 
focus our attention on the n = 2 case where each of the TV 
bosonic atoms has only two available states, denoted as j" and I 
respectively: in this case, the symmetric border rank, the sym- 
metric rank and the decomposition can be obtained using the 
algorithm reviewed in Sec l3.2l 



4.1 Spin Squeezed States 

In the last years, spin squeezed states (SSS) have attracted quite 
some attention as a tool to improve the sensitivity of inter- 
ferometric device for precision measurements, as they allow 
to overcome the so-called standard quantum limit imposed by 
shot noise M27U28II . Experimental demonstrations using differ- 
ent physical systems have been reported in |30 l 3 1U321 . 

To generate a spin squeezed state, one can start from a co- 
herent spin state (CSS), i.e. a symmetric rank-1 Hartree state 
where all atoms share the same wavefunction. In the present 
case where the single atom Hilbert space is bidimensional, Har- 
tree states can be parametrized in terms of a macroscopic spin 
vector moving on the surface of the so-called Bloch sphere. Its 
direction is specified by two angles 8, if with < 6 < ir and 
< Lp < 2?r, 



|coh : 6,ip) = 



/TV! 



cos(0/2) a\ + e lv sin(0/2) at 



N 



vac) 



(23) 
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the 9 — (6 = it) case with the macroscopic spin pointing 
along (opposite to) the quantization axis z corresponds to all 
atoms being in the j" (X) state. A macroscopic spin along the 
"horizontal" plane orthogonal to the quantization axis corre- 
sponds to equally populated j" and \, states; the phase ip corre- 
sponds to the angle made by the macroscopic spin and the x 



As proposed in [28], SSS can be obtained from a 6 = ir/2, 
Lp = CSS by means of a one-axis twisting of this state via the 
self -phase modulation Hamiltonian, 

\SSS : fx) = exp(-i(xS*/2) |coh : 6 = n/2,ip = 0). (24) 

In elementary quantum mechanics, the total spin operator S z is 

(i) 

defined as the sum of the corresponding Si operators acting 
on the z-th atom and leaving the other atoms unperturbed; the f 

and I states of the i-th atom are eigenstates of Si with eigen- 
value ±1/2. In a second quantization formalism, the total spin 
operator S z has the simple form 



W) 



(25) 



The first step in view of computing the rank of the | SSS) 
state is to write it in a second quantization form as a polyno- 
mial in the creation operators a j ^. Expanding the coherent spin 
state ( 1231 and then applying the one-axis twist operators as in 
d24l) . one easily obtains the form 



In order for this equality to hold, one has to have fi(jg —j a ) = 
(mod 2tt) for at least one pair j a ^ jp. 

As discussed in Sec 13. 21 the state has symmetric border 
rank equal to r if and only if the rank of the r-th catalecticant 
matrix CN- r ,r(SSS) is not maximal (i.e. smaller than r + 1) 
and the rank of all the previous r'-th catalecticant matrices 
Cjsr- r ',r' (SSS) with r' < r was maximal (i.e. equal to r' + 1). 
Assuming that this latter condition is satisfied, the symmetric 
border rank is r if for all possible choices for ji, . . . ,j r +u 
there exists at least one pair j a ^ jp such that 



nUp-Ja) = 0(mod27r). 



(30) 



[2^1 corres P on ds to 



If \x is an irrational multiple of 2w, this condition can never 
be satisfied, so the symmetric border rank of the spin squeezed 
state has the generic value qn — \^^- 
have maximal symmetric border rank. 

If /j, is a rational multiple of 2ir, i.e. \i = 2np/q (p, q are 
coprime integers) and q < r, the pigeon hole principle guar- 
antees that for any choice of row indices < ji, ■ . ■ ,j r +i < 
(N — r + 1), there always exists a pair j a 7^ jp such that 
n(j a — jp) = (mod27r). On the other hand, for r < q there 
exists a set j\ — 1, . . . ,j r +i = r + 1 such that all pairs 
fi(j a — ip) ^ 0(mod27r). Combining these two arguments 
into the algorithm of Sec l3.2l we conclude that the symmetric 
border rank of the spin squeezed state with squeezing param- 
eter 11 = 2irp/q is equal to q if q < and equal to if 

q > 9n- 



\SSS : fx) 

N 



vww. 

(^) e***- N t*f'* (a t ) N ~ k (a,)" |vac). (26) 



fc=0 



The key element of the algorithm discussed in Sec l3.2l are the 
catalecticant matrices CN- r ,r(SSS) associated to the polyno- 
mial: in our case, their entries for l<fc<A^ — r + 1 and 
1 < j < r + 1 are equal to 



1 



Cfc, 



y/2 N Nl 



exp 



N 



(27) 



To determine whether the SSS has symmetric border rank r 
one has to test the r + 1-minors of the r-th catalecticant matrix 

C]y— r , r (SSS). 

Let < ji, . . . ,j r +i < (N — r + 1) be the row indices 
of the rows that are involved in a given minor. Straightforward 
algebraic manipulation shows that the minor is zero if and only 
if the Vandermonde determinant l39ll 



K 

ZJ 6'- 

O 

pa 



1/3 



(2 



(X/27U 

Fig. 1. Numerical plot of the symmetric rank and symmetric bor- 
der rank of spin squeezed states of the form d24l l as a function of the 
squeezing parameter (x. Total number of particles N = 13. 



I e -wi . . . g-'Mir 
I e ->Hi . . . g-W2r 



I g-Wr+1 . . . g-Wr+ir 

vanishes, i.e. if and only if 



n 

l<j a <3><(r+l) 



(28) 



■)=0. 



(29) 



In order to determine the value of the symmetric rank, we 
have to consider the polynomial q(to,tx) defined in (fT~5b that 
is associated to vectors in the null space of CN-r.r(SSS) for 
r = sbr(SSS) 

For [x = 2-Kp/q and q < g^, direct inspection shows that 
the (r + 1) -component (r = q) vector 



(l,0,...,0,-e OTp<? e- OTpAr )' 



(31) 
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always belongs to the null space of the r-th catalecticant matrix 
CN-r.r(SSS), so that the r roots of the associated polynomial 



(32) 



are always distinct. This shows that the symmetric rank is al- 
ways equal to the symmetric border rank: 



sr{SSS) = sbr(SSS). 



(33) 



The situation is more delicate in the general case when [i 
is an irrational multiple of 27r or [i = 2irp/q with q > g^, so 
that the symmetric border rank has the maximal, generic value 
sbr = <7jv. If N is even, the only possibility for the symmet- 
ric rank is to be sr = sbr = g^- If N is odd, one has to 
assess whether sr = g^ (the generic case) or sr = gx + 1. 
To this purpose, one has to repeat the procedure of calcula- 
tion the roots of the polynomial associated with the null space 
of the r = <7jv + 1 catalecticant matrix. By means of a MAT- 
LAB/OCTAVE code implementing the Algorithm described in 
Sec l3.2l we have checked that the roots are always distinct, so 
that sr = sbr = g^. These results on dependence of the sym- 
metric rank of SSS on the squeezing parameter [i are summa- 
rized in FigQ] 

Let us now move to the decomposition. For /i = 2-7r/g, 
the roots of the polynomial d32l are complex numbers of mag- 
nitude 1 and equidistant phases by 2ir/q: this remarkable fact 
recovers a result originally found in 11331 that the spin squeezed 
state with /i = 2up/q can be written as a Schrodinger-cat- 
like linear superposition of q macroscopically different coher- 
ent spin states of the form (1231 with 6 = n/2 and equidistant 
phases (j> = </>o + 2ixj /q, j = 0, . . . , q — 1; the initial phase 
is 0o = unless both p and q + N are odd, in which case is 
<Po = K/q. 

The direction of the macroscopic spin is therefore always 
oriented along the plane orthogonal to the quantization axis and 
makes equidistant angles <j> with the x direction. A numerical 
solution of the linear system at the step 6 of the algorithm pre- 
sented in Sec 13. 21 shows that the weigths A; of the q coherent 
spin states have equal magnitude. These facts are summarized 
by the remarkable writing of the spin squeezed state in the form 

\SSS : fjL = 2irp/q) = 

= V — |coh : 9 = n/2, <p = O + 2nj/q) (34) 
H ^ 



with suitably chosen phases Q'2 [33 



4.2 Fock States 

In the previous subsection we have seen that the symmetric 
rank and symmetric border rank are equal for spin squeezed 
states for all values of p. Here we illustrate a class of physically 
relevant states for which the symmetric rank and symmetric 
border rank are different. Having this distinction clear in mind 
is essential if the concept of Schmidt rank is to be used as a 
measure of entanglement l22l . 



In the mathematical literature, it is a well known fact that 
the symmetric rank and the symmetric border rank of bivari- 
ate monomials usually have different values: straightforward 
inspection of the form of the catalecticant matrices and of the 
associated polynomial (JT3J shows that the symmetric rank of 
a monomial of the form x k y , with k < N/2, is equal 
to sbr = k + 1, while the symmetric rank is equal to sr = 
N - k + 1 > sbr. 

The physical states corresponding to monomials are Fock 
states (FS) 



IV') = (at) 



t(JV-fe) 



(%) 



vac) 



(35) 



where exactly k particles sit in the I state and the remaining 
N — k particles sit in the j" state. The results on the symmet- 
ric rank and the symmetric border rank of monomials quoted 
above imply that a Fock state can be decomposed either as the 
exact sum of sr = N — k + 1 CSS states or as the limit of the 
sum of a smaller (equal only if k = N/2 and N even) number 
sbr = k + 1 < sr of CSS states. 

The decomposition as an exact sum of sr CSS states has 
the form 



(a t )t^(a i )t fc |vac) = 



1 



(N 



1)( 



N-k 
3=0 



N 



(36) 



where the relative phases of the f, J, amplitudes in the sr CSS 
states have been carefully chosen in a way to have a complete 

destructive interference on all monomials a^~ k a| fc with k' ^ 
k. 

The decomposition as the e — > limit of a sum of sbr CSS 
states has a similar structure: 



(4) w - fc (4) fc |vac) 



1 



x lim —r 

e->0 e k 



E(4 + ^-M 



3=0 



N 



k a. 



I vac). (37) 



In contrast to Q61) . only the monomials a| fc with k' < k vanish 
by destructive interference, while the ones with k' > k disap- 
pear when the 640 limit is taken. 



5 Conclusions and perspectives 

In this paper we have reported an example of application of al- 
gebraic geometry tools to a simple model of quantum physics. 
We have calculated the symmetric rank and the decomposi- 
tion of spin squeezed states and Fock states of an assembly of 
N indistinguishable two-level bosonic atoms: the rank of spin 
squeezed states shows interesting features when the squeezing 
parameter has a rational value (in suitable units) so that the 
state can be decomposed as the superposition of a few macro- 
scopically different coherent spin states. Fock states have the 
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peculiar property that their symmetric rank and symmetric bor- 
der rank have different values. 

A most challenging future perspective is to generalize the 
decomposition algorithms to the case of a large number N of 
atoms and a large number m of available states and then com- 
bine them with a numerical simulation of the time-evolution 
of an interacting system. First attempts in this direction used 
the semiclassical representations of a quantum field l42l or 
a reformulation of the A^-body problem in terms of stochas- 
tic Hartree l43l (Hartree-Fock R4l ) wavefunctions for systems 
of bosonic (fermionic) particles, but in many practical appli- 
cations to strongly correlated many body systems this method 
lead to numerically untractable problems. 

From a mathematical point of view, the computation of the 
symmetric border rank of a homogeneous polynomial of degree 
N in n > 2 variables (i.e. A" bosonic atoms in n > 2 avail- 
able states) is a still open problem as the equations defining the 
corresponding varieties are not yet understood (see e.g. 11401 
and references therein). Nevertheless, algorithms to compute 
the symmetric rank and the decomposition have been recently 
developed (26] using a generalized version of the catalecticant 
matrices, the so-called Henkel matrices. Their application to 
quantum physics problems is a natural next step for our re- 
search. 

The situation for the skew-symmetric tensors describing 
systems of N fermions is somehow specular: a criterion to 
determine the skew-symmetric border rank of a given skew- 
symmetric tensor of any N was recently developed 1411 . On 
the other hand, we are not aware of any algorithm to compute 
the skew-symmetric rank and the actual decomposition except 
for the first meaningful case A" = 2 where the tensor is a ma- 
trix. 
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